Simplest model to study reentrance in physical systems 
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We numerically investigate the necessary ingredients for reentrant behavior in the phase diagram of physical 
systems. Studies on the possibly simplest model that exhibits reentrance, the two-dimensional random bond 
Ising model, show that reentrant behavior is generic whenever frustration is present in the model. For both 
discrete and continuous disorder distributions, the phase diagram in the disorder-temperature plane is found to 
be reentrant, where for some disorder strengths a paramagnetic phase exists at both high and low temperatures, 
but an ordered ferromagnetic phase exists for intermediate temperatures. 



PACS numbers: 75.10.Hk, 05.70.Fh, 64.60.Cn, 64.60.an 

Reentrance in a thermodynamic phase diagram presents a 
counterintuitive scenario: one phase exists inside some closed 
temperature range, with transitions to the same second phase 
at both higher and lower temperatures. If one of these phases 
is more ordered than the other one, then one of the phase 
transitions must violate the intuition that lower temperature 
phases are more ordered than higher temperature phases. In 
the context of solid-liquid transitions, this is known as "in- 
verse melting" or "inverse freezing." A variety of materi- 
als have been shown to have reentrant phase diagrams. For 
example, Rochelle salt was found to be a ferroelectric with 
two Curie points; the ordered ferroelectric phase occurs only 
between these two temperatures More recently, similar 
phase diagrams have been seen for superconducting vortices 
El O, liquid crystals m, miscibility in solutions fS), poly- 
meric materials [61, ferromagnetism in semiconductors Q, 
denaturation of DNA |8|, and many other systems iP.TOl. 

Theoretically, a number of model systems have been found 
with reentrant phase diagrams. The Ising model on a Union 
Jack lattice with frustrated anisotropic interactions was shown 
to be reentrant in a narrow region of parameter space ifTTlfTZl . 
In the fully-frustrated Villain model ||13| the ground state of 
the system is seen to be disordered, while the low-lying ex- 
cited states favor ferromagnetic ordering; in this case the or- 
dering is due to the emergence of ferrimagnetism 1 14). Reen- 
trant ferromagnetism has also been seen in models of semi- 
conductors because the carrier density increases with temper- 
ature 1 15 1; in models of solid hydrogen, reentrance is due to 
quantum fluctuations 1 16|. Most prominently, frustrated spin- 
glass models have been shown to include the complexity nec- 
essary to describe rather generic reentrant scenarios ifTTHSTl . 
albeit for models with complex Hamiltonians. 

Here we show numerically that the two-dimensional 
random-bond Ising model (RBIM) — possibly the simplest 
disordered spin model with frustration — generically possesses 
reentrance in its phase diagram for both discrete and continu- 
ous disorder distributions. It is given by the Hamiltonian 

with an L X L square toroidal grid 1221 of Ising spins {s,; e 
±1} and quenched nearest-neighbor random couplings 



These random couplings are most commonly chosen from ei- 
ther a bimodal (± J) or Gaussian distribution. We emphasize 
that this Hamiltonian is much simpler than the disordered spin 
models typically used to study inverse freezing (e.g., the "sim- 
plest" model of inverse freezing Iil9il ). Furthermore, the phase 
diagram of this model is uncomplicated. At any finite temper- 
ature, with a disorder distribution biased toward ferromagnetic 
interactions, only ferromagnetic and paramagnetic phases are 
possible 1231 . Because reentrance scenarios often occur in 
complex phase diagrams with many distinct phases, such a 
clean example is difficult to find in reentrant models. 

The random-bond Ising model. — The RBIM is widely 
studied as a standard model of disordered systems; it is one 
of the simplest models which includes the disorder and frus- 
tration needed to exhibit a complex glassy behavior at low 
temperatures [24] . These same elements, especially frustra- 
tion, are necessary for the reentrant behavior in the phase 
diagram seen here. Furthermore, the model and cousins on 
more complex lattice geometries are of paramount importance 
across disciplines and more recently have found widespread 
use in the computation of the error stability of topologically- 
protected quantum computing proposals l25lj27l . For quan- 
tum computing applications, the order-disorder transition cor- 
responds to the maximum error rate at which quantum oper- 
ations may be performed with high fidelity. This relation is 
particularly apparent because the quantity we use to identify 
the phase transition in the RBIM is closely related to the error 
rate in the quantum computing proposals. 

We study the phase diagram for this model with both dis- 
crete and continuous bond disorder { J^ }. For the discrete ± J 
distribution, bond values are chosen according to 

P( J,,) = p5{J,j - J) + (1 - p)5{.J,j + J). (2) 

The pure Ising model is recovered forp 1, while the strong- 
disorder case is typically studied for p = 0.5. The disorder 
strength q ^ 1 — p gives the deviation from a pure ferromag- 
net, but note that the free energy F of the model is symmet- 
ric under reflection about p — 0.5, so F{T,p) = F{T,q). 
Nishimori has shown that, due to extra symmetries of the 
problem, a number of quantities, such as the internal energy 
of the system, may be computed exactly when the equality 
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(1 -p)/p = exp(-2 J/T) holds f2^. This equahty is called 
the Nishimori line. The location of the phase transition on 
the Nishimori line, (T*,p*), has been identified with a mul- 
ticritical point Of interest here is that the magnetization 
for a given value of p must be greatest on the Nishimori line, 
so that the ferromagnetic phase must not exist at any temper- 
ature for p < p* |28|. This implies that the phase diagram 
below the Nishimori line is either vertical or reentrant. It has 
been argued analytically that the vertical case is the correct 
one |30, 31 1, although numerical studies with ±J disorder 
suggest a reentrant phase diagram 1.32 - 1-35 J . 

For the Gaussian distribution, bonds are chosen from 
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exp 



(3) 



The pure Ising model is recovered for J ^ 0, while the strong 
disorder case occurs when J = 1 and J — 0. The free energy 
of the model is symmetric under reflection about J = 0. It is 
customary 13611371 to define a disorder strength parameter r = 
J/ J. In the case of Gaussian disorder, the Nishimori line 1281 
is given by J^/T = J. Here, the multicritical point (T*, r*) 
is also the largest value of r for which a ferromagnetic phase 
may exist, so that the phase diagram for Gaussian disorder is 
also expected to be either vertical or reentrant. 

Unfortunately, the exact location of the multicritical point 
has not been calculated for either disorder distribution; nu- 
merical estimates for ± J disorder have given values of q* = 
0.1094(2) 1381, 0.1093(2) m, and 0.10917(3) JH, while 
estimates for Gaussian disorder include r* = 0.97945(4) 
BOl and 0.9811(3), where the latter number was quoted as 
1/r* = 1.0193(3) im. Exact efficient ground state algo- 
rithms exist for studying this model at zero temperature |42), 
although the use of such techniques is complicated because 
they typically do not handle degeneracy well. Nevertheless, 
the values obtained are in the vicinity of qc{T = 0) = 0.103 
lEIllMl and rc(T = 0) = 0.970(2) [37|, and clearly are not 
consistent with the values at the multicritical point. 

One could still imagine a scenario where the T = behav- 
ior is significantly different from all finite-temperature behav- 
iors: this is the case for the strong-disorder region of the phase 
diagram, where the T — spin-glass "phase" does not exist at 
nonzero temperatures ll23l . Low -but-finite-temperature simu- 
lations are typically quite difficult due to long equilibration 
times, and the difference between the numerical estimates is 
quite small, so few points in between have been probed. The 
numerical measurement of the location of the phase transition 
at finite temperatures below the Nishimori line 1351 l39l has 
required thorough finite-size scaling extrapolation. Neverthe- 
less, computations of Pc{T) for two points below the Nishi- 
mori line have shown intermediate results I.35J . Here, our 
more efficient technique allows us to see the reentrance far 
more clearly than in previous studies. 

Simulation details. — To investigate the disorder- 
temperature phase diagram, we numerically identify the 
parameters that give the ferromagnet-paramagnet phase 
transition using a Pfaffian technique. Because this phase 



transition is of second order, the magnetization, which is the 
natural order parameter, is continuous. Quantities such as 
the Binder ratio |43 | have been developed to most precisely 
determine the location of a phase transition. Here, we use a 
different quantity, which is more easily computed when the 
partition function is directly available. 

We start be defining an extended Hamiltonian (see 
Ref |44|) where the boundary conditions are allowed to vary, 
with both periodic and antiperiodic cases being allowed. The 
extended Hamiltonian is given hy H* = — Jij^iSjaij, 
with (Tij = 1 except on one vertical column of horizontal 
bonds where aij ~ a-y and one horizontal row of vertical 
bonds, where aij — ah', the ah^y = ±1. Defining a config- 
uration {{s,;}, CTy, ail} of the system now requires specifying 
the spin values as well as the boundary conditions. The parti- 
tion function, Z — J2{s =±i} exp(— H/T), may also be ex- 
tended to Z* = E{.,=±i}..„=±i,.,=±i cxp(-HVr). Us- 
ing a Pfaffian technique, Z* may be directly evaluated from 
the same computation that produces Z with no additional 
computational effort. In this extended system we can compute 
the probability »c that the boundary conditions are periodic in 
both directions, given by 



Z/Z*. 



(4) 



When the boundary conditions in a direction are changed from 
periodic to antiperiodic a system-spanning domain wall of 
length scale L is imposed across the system. In the ferro- 
magnetic phase, a system-spanning domain wall is energeti- 
cally unfavorable: the periodic -periodic case, ay = a^ = 1, 
has lower free energy by an amount proportional to L and 
therefore sc = 1 as i — > oo. In the paramagnetic phase, the 
boundary conditions do not significantly influence the thermo- 
dynamics of the system, and all four cases have equal weight, 
so that >K = 1/4 as L — oo. At the phase transition, >k ap- 
proaches a constant independent of L (up to finite-size correc- 
tions), i.e., »c ~' ZiL^/^iT-Tc)]. This measure is somewhat 
more sensitive than the free energy of a system-spanning do- 
main wall because it allows for system-spanning domain walls 
in either or both directions. Note also that 1 — sc is equiva- 
lent to the error rate discussed in Ref. [33 J. Thus the method 
we introduce here can be applied generically to the study of 
topologically-protected quantum computing proposals. 

We have computed sc to investigate the phase diagram of 
the random-bond Ising model using a Pfaffian technique. With 
periodic boundary conditions, one may exactly compute Z by 
summing the Pfaffians of four related Kasteleyn matrices |45 1. 
The partition function of a system of N = Lx L spins may be 
computed in 0{N^^'^) operations, allowing for exact partition 
function evaluation in systems up to L = 512. By choosing 
different signs for the terms in the sum, it is also possible, 
without computing any additional Pfaffians, to compute Z for 
all four boundary conditions [451 . Thus by computing Z, we 
obtain Z*, and therefore jk for no extra effort. 

Results. — Figure [T] shows the order probability >k as a 
function of temperature T and disorder strength q for the 
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FIG. 1: (Color online) Order probability function of tempera- 

ture T and disorder strength q for the RBIM with ± J interactions for 
system sizes L = 16, 32, and 64 (shallowest to steepest). Finite-size 
scaling corrections are small, i.e., the surfaces cross cleanly at the 
phase boundary (shown in the projection onto the plane). Error bars 
are smaller than the thickness of the surfaces. 



random-bond Ising model with ± J interactions and different 
system sizes. The data do cross at a line (see projection) that 
corresponds to the phase boundary, thus illustrating that the 
approach used works well. To extract the best estimate of the 
critical temperature Tc for a given value of q we perform a 
finite-size scaling of the data with and ly as free parameters. 
After performing a Levenberg-Marquard minimization of the 
chi^ of the best fit to a third-order polynomial we estimate sta- 
tistical errorbars by wrapping the process in a bootstrap anal- 
ysis. The phase diagram for ±J disorder is shown in Fig. |2] 
and clearly shows reentrant behavior. We have these scaling 
collapses for data at fixed q (squares) and T (circles). To fur- 
ther highlight the reentrant behavior, in Fig.|3]we show the or- 
der probability as a function of temperature T for q ~ 0.107 
(vertical line in Fig. |2]i. The data show two crossings, there- 
fore clearly indicating that the phase diagram is paramagnet- 
ferromagnet-paramagnet with two distinct transitions. 

Finally, in Figs.|4]and[5]we show data for Gaussian disorder. 
Reentrance is clearly present, albeit much weaker than for the 
±J case: the ratio rc{T = 0)/r* « 0.99 is much closer to 
1 than qc{T = 0)/q* « 0.94 for the ±J case. These results 
show clearly that reentrance is a generic feature of this model 
when disorder and frustration are present. The fact that the 
case with Gaussian disorder has a much weaker effect sug- 
gests that the ground-state entropy might play a role but is 
not strictly necessary. In fact, studying a model where a con- 
tinuous transition between the ± J and Gaussian cases can be 
tuned 1 47 1 might help in elucidating this behavior, but it would 
be computationally extremely expensive. This tuning could 
change the magnitude of the reentrance, but it appears that the 
phase transitions are in the same universality class: the criti- 
cal exponent v is consistent with z/ « 1.5 for all points below 
the Nishimori line. For ± J disorder, we find an aggregate 
1/ = 1.49(4) while for Gaussian disorder, = 1.52(5), in line 
with previous studies Il33ti35l [37l 139^411 . 

Summary and Discussion. — We have shown that the 
random-bond Ising model in two dimensions — a simply posed 
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FIG. 2: (Color online) Phase diagram of the two-dimensional 
random-bond Ising model with ±J interactions [Eq. |2|]. The 
shaded region is ferromagnetic, while the white region is paramag- 
netic. The boundary shown between these two regions is a guide 
to the eye. Some phase boundary points are computed by doing a 
scaling collapse varying T (squares), while others are from a scaling 
collapse varying q (circles). Statistical error bars which are smaller 
than the symbol size are not visible. For q < 0.05, we use the results 
from Ref. |46|, which used a related technique. The black circle for 
g = is the exact result for the pure Ising model. The point on 
the Nishimori line is from Ref. [35 1 and the T — point where the 
boundary touches the axis is from Ref. |34). 
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FIG. 3: (Color online) Order probability function of tempera- 

ture for q — 0. 107 (vertical line in Fig.|2]l and ± J disorder. The data 
show two crossings, illustrating the existence of two transitions and 
therefore reentrance in the phase diagram. 



model with only two phases — generically possesses reen- 
trance in its phase diagram. The disorder and frustration 
present in this model are responsible for this counterintuitive 
result. This is likely related to the "order by disorder" seen in 
the Villain fully frustrated model 114|. In the fully-frustrated 
case, ferromagnetic strips are completely decoupled from one 
another in the ground state, but the low-lying excitations have 
a weak ferromagnetic interaction among the strips. In the 
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FIG. 4: (Color online) Phase diagram of the two-dimensional 
random-bond Ising model with Gaussian disorder [Eq. l|3j]. The 
shaded region is ferromagnetic, while the white region is paramag- 
netic. The black points on the Nishimori line are from Refs. 140J 
(left) and |41 1 (right), and the point at T = is from Ref. IST). The 
boundary shown between these two regions is a guide to the eye. 




FIG. 5: (Color online) Order probability >k as a function of temper- 
ature for r = 0.976 (vertical line in Fig.|4]( and Gaussian disorder. 
Although the approach to the thermodynamic limit is slow, the data 
clearly show two crossings, illustrating the existence of two transi- 
tions and therefore reentrance in the phase diagram. 



RBIM, the ground state consists of ferromagnetic domains of 
a size and energy scale set by the disorder strength. When 
these are dense enough to percolate throughout the system, 
ferromagnetic ordering will cease at T = 0. However, for a 
range of parameters these domains can be coupled strongly 
enough in the low-lying excitations to produce the reentrant 
behavior in this model. It should be noted that the reentrance 
scenario shown here is particular to two space dimensions; 
the random bond Ising model in higher dimensions has a low- 
temperature spin glass phase so that if the ferromagnetic phase 
only exists for intermediate temperatures, the low-temperature 
phase would be a spin-glass phase and not the same paramag- 
netic phase found at high temperatures. 
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